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We present a local probe study of the magnetic superconductor, ErNi 2 B 2 C, using magnetic force 
microscopy at sub-Kelvin temperatures. ErNi 2 B 2 C is an ideal system to explore the effects of con¬ 
comitant superconductivity and ferromagnetism. At 500 mK, far below the transition to a weakly 
ferromagnetic state, we directly observe a structured magnetic background on the micrometer scale. 

We determine spatially resolved absolute values of the magnetic penetration depth A and study its 
temperature dependence as the system undergoes magnetic phase transitions from paramagnetic 


to antiferromagnetic, and to weak ferromagnetic, 
dition, we estimate the absolute pinning force of 
temperature dependence as well, and discuss the 
formation. 

PACS numbers: 74.25.Wx, 74.25.Ha, 68.37.Rt 

I. INTRODUCTION 

The interplay between magnetism and superconduc¬ 
tivity has caught the attention of many condensed mat¬ 
ter scientists. Once thought to be mutually exclusive, 
several materials have since emerged in which a coex¬ 
istence of magnetic order and superconductivity can be 
found.®® This is contrasted, e.g., by high temperature 
cuprate superconductors; although a low temperature 
rare-earth-related magnetic ordering may occur within 
the superconducting phase,® there is a clear separation 
between the Cu-based strong antiferromagnetic and the 
superconducting phases.® Nevertheless, spin fluctuations 
in the Cu-0 planes that arise from the melting of the an- 
tiferromagnetically ordered phase upon doping seem to 
support and mediate superconductivity. 5 6 

The possibility of studying materials that exhibit mag¬ 
netic order and superconductivity simultaneously could 
shed a new light on the importance of magnetic correla¬ 
tions for the pairing mechanism of Cooper pairs in uncon¬ 
ventional superconductors. In particular, our insight into 
high temperature superconductivity could greatly bene¬ 
fit from a detailed study of the interplay among intrinsic 
magnetism, penetration depth A, and coherence length 
£■ 

From an applied point of view, the pinning of 
Abrikosov vortices in type-II superconductors is a cen¬ 
tral issue, as increased pinning forces Fp can greatly 
enhance both critical currents J c and critical magnetic 
fields H c2 .® The presence of an intrinsic magnetic field 
(with Bi n t < H c 2 ) can have a strong effect on Fp and 
may thus present an alternative to artificially introduced 
pinning centers. 


, all within the superconducting regime. In ad- 
Abrikosov vortices, which shows a position- and 
possibility of the purported spontaneous vortex 


The rare earth borocarbides i7Ni 2 B 2 C (with R = Dy, 
Ho, Er, and Tm) exhibit magnetism within the super¬ 
conducting phase. 1 Conspicuously, their superconducting 
and magnetic transition temperatures are comparable, 
suggesting similar energy scales for magnetism and super¬ 
conductivity. ^ The observation of a pronounced isotope 
effect in these intermetallic compounds points towards 
conventional phonon mediated superconductivity.®^ Of 
these borocarbides, ErNi 2 B 2 C has a relatively high su¬ 
perconducting transition temperature of T c = 10.5 K 
and develops antiferromagnetic order around T/v = 6 
K.®^ Remarkably, below Twfm = 2.3 K a weak fer¬ 
romagnetic phase can be detected with a net magnetic 
moment of 0.39/ie/Er and a resultin g intr insic magnetic 
field of about B int = 500 — 700 G ? 11 * 12 ! Previous tun¬ 
nel diode oscillator (TDO) experiments focusing on the 
magnetic penetration depth A uncovered a pronounced 
temperature dependence with clear features at T/v and 
Twfm , therefore indicating a strong impact of the mag¬ 
netic properties on A, as well as fingerprints of the for¬ 
mation of a spontaneous vortex lattice (i.e. induced 
by intrinsic magnetism) below Twfm f 13 At the same 
time, Bitter decoration experiments 14 and scanning Hall 
probe measurement^! revealed a microscopic variation 
of the intrinsic magnetic field emerging below T/v in the 
form of well-ordered, micrometer sized magnetic stripes, 
along which vortices tend to order. A spontaneous vor¬ 
tex phase in the weak ferromagnetic regime was not ob¬ 
served, though. However, the Bitter decoration exper¬ 
iments were conducted down to 1.9 K, which is only 
slightly below Twfm , and the spatial resolution of scan¬ 
ning Hall probes is rather limited. In fact, until now no 
clear consensus on the issue of spontaneous vortex lat- 
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tice formation in ErNi 2 B 2 C has been reached, partly due 
to its low ferromagnetic transition temperature. This 
open issue together with the microscopic variation of 
the magnetic field calls for a local probe investigation 
of ErNi 2 B 2 C at sub-Kelvin temperatures. 

Here we present a magnetic force microscopy (MFM) 
study on ErNi 2 B 2 C at He-3 temperatures and with a 
high spatial resolution. We observe magnetic stripe¬ 
like features at T = 500 mK, that arise from twin do¬ 
main boundaries. Their pronounced temperature behav¬ 
ior supports the spontaneous vortex formation scenario. 
Using the magnetic moment of the tip we estimate the 
pinning force of a single vortex. In addition, a spatially 
resolved determination of the penetration depth A is pre¬ 
sented, which reveals a close relation between A and the 
microscopically varying intrinsic magnetism. 


II. EXPERIMENTAL DETAILS 

A single crysta l of E rNi 2 B 2 C was grown via Ni 2 B flux 
growth technique. 10 15 The sample’s magnetic properties 
have been measured to determine T c using an MPMS 
superconducting quantum interference device (Quantum 
Design, Inc.). The crystal orientation was determined by 
polarized Raman scattering experiments, using the ori¬ 
entation dependence of the B\ g phonon mode intensity.® 
Prior to MFM measurements, the sample was mechan¬ 
ically polished to obtain a fresh surface within the ab 
plane with a roughness as small as 5 nm. MFM mea¬ 
surements were performed with a unique home-built He- 
3 MFM probe with a vector magnetic field capability (2T 
- 2T - 9T along the x-y-z direction) and a temperature 
range of 320 mK - 300 K. 17 All experiments were carried 
out with no magnetic fields applied and with a commer¬ 
cially available tip. 18 

III. EXPERIMENTAL RESULTS & DISCUSSION 

We first discuss the origin of magnetic stripes and the 
local variation of the magnetic penetration depth, be¬ 
fore moving on to the pinning force. Synchrotron x-ray 
scattering experiments on ErNi 2 B 2 C revealed a struc¬ 
tural tetragonal-to-orthorhombic distortion upon enter¬ 
ing the antiferromagnetic phase at T/v = 6 K. This leads 
to the formation of twin domain boundaries below T/v 
along [110] and [110] p- 9 which are already inhabited by 
a ferromagnetic component aligned along the c axis.® 
In small applied magnetic fields (of the order of 10 G) 
vortices were found along the boundaries, thus forming 
stripes with roughly 5 ym spacing. 21 The transition at 
Twfm is driven by magneto-elastic coupling, which en¬ 
hances Ruderman-Kittel-Kasuya-Yosida (RKKY) inter¬ 
actions among the Er 4/ electrons, and thereby a long- 
range ordered weak ferromagnetism is induced in the 
bulk, as evidenced by inelastic neutron scattering ex- 
perimentsP 2 In our MFM data [see Fig. 1(a)] a struc¬ 


tured magnetic background can be clearly seen in the 
form of bright stripes with a separation of « 7.5 fi m 
at T = 500 mK, i.e., deep within both the supercon¬ 
ducting and the weak ferromagnetic phase. The stripes 
run from the lower right corner to the upper left corner 
[highlighted by dashed lines in Fig. 1(a)] and at 45° from 
the crystallographic a and b axes. We therefore ascribe 
the stripes to regions of increased ferromagnetism. Their 
orientation along the [110] or [110] directi on su ggests a 
coincidence with twin domain boundaries P - I-l l Indeed, 
magnetic stripes of the same orientation and comparable 
spacing have been observed in Bitter decoration experi¬ 
ments for T < T/v and attributed to twin boundaries. ®- 20 - 1 
However, the magnetic stripes in our experiment evi¬ 
dence a considerably different temperature behavior, as 
discussed further below. 

In order to probe the influence of the varying internal 
magnetic field on the superconductivity in the sample, 
we performed local measurements of A across two stripes 
along the blue arrow in Fig. 1(a). The determination of 
A can be done by slowly bringing the tip to the surface at 
a fixed sample position and measuring the Meissner force 
gradient acting upon the magnetic tip. By comparing the 
Meissner force curves to those of a well-known standard 
sample (a Nb film), one can obtain the absolute values 
of A. 23 In Fig. 1(b) we plot two Meissner force curves 
measured at the beginning (dashed black) and at the 
end (solid green) of the arrow. A clear shift is observed, 
indicating a spatial variation in A. In the inset of Fig. 
1(b) several absolute A values at different positions (open 
squares) are plotted. These values, measured under the 
same conditions, vary by more than 20% within a few 
micrometers, while following closely the changes in the 
intrinsic ferromagnetism due to the stripe formation.® 
This is to be expected, as the magnetic penetration depth 
will be enhanced in the presence of an increased intrinsic 
magnetic field. In fact, we can reinterpret our MFM 
images in the superconducting state as maps of changes 
in A. 

The temperature dependence of A [see Fig. 1(c)] can 
be well-described by a BCS-type behavior, i.e., A(T) = 

yji X (T/T y ’ w ^ ere ^(0) penetration depth at T = 0 

K. The fit yields A(0) « 190 nm and T c « 9.1 K. 
Our directly obtained penetration depth A(0) is larger 
than those previously reported from bulk probe experi¬ 
ments P^ 25 while T c is slightly below the reported value of 
10.5 K. Aside from the different experimental approaches, 
it is known that the T c of as-grown samples is lower than 
that of annealed samples.® Indeed, bulk magnetization 
measurements of our as-grown sample yield a T c of 9.9 K 
[see red line in Fig. 1(c)]. Another important consider¬ 
ation is a possible deviation of ErNi 2 B 2 C from the BCS 
behavior. In fact, a pronounced non-BCS-type behav¬ 
ior is apparent below T c in the highly resolved tempera¬ 
ture dependence of the TDO study,® which is related to 
the different magnetic phases within the superconduct¬ 
ing regime. Since our temperature dependence is more 
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FIG. 1. (Color online) (a) A 15 fine i x 15 /im MFM image obtained at T — 500 mK and 300 nm tip-sample distance. The blue 
arrow marks a line of 6.5 /im length along which spatially resolved A measurements were performed. The dashed lines outline 
the magnetic stripes, (b) Meissner force curves obtained at the beginning [i.e. black dot in (a)] and at the end of the arrow 
[i.e. green dot in (a)]. The inset shows local A values along the blue arrow in (a), (c) Temperature dependence of A obtained 
at a single position (blue spheres) and the magnetic moment measured at H = 30 Oe (solid red line). The dashed blue line is 
a fit according to the BCS theory. 


coarse, a deviation from the BCS behavior could be ob¬ 
scured. 

In addition to the magnetic stripes, we clearly see four 
separate vortices in Fig. 1(a). They are aligned antipar¬ 
allel to the magnetic stripes, as indicated by the MFM 
contrast. Furthermore, they remain above Twfm and 
even above T/v (as discussed below). Therefore, they 
exist independently of any intrinsic magnetic property. 
We ascribe these features to vortices induced by a weak 
external stray field. As this stray field is present regard¬ 
less of the temperature (i.e., already above T c , in con¬ 
trast to the weak intrinsic ferromagnetism), its flux lines 
will be trapped inside ErNi 2 B 2 C upon cooling through 
T c . The observation of n = 4 vortices on a scan area 
of A = 15 x 15 /im 2 could already be accounted for by 
Hstray = ~ 0.37 Oe, considering that each vortex 

corresponds to a flux quantum of $o. Although this esti¬ 
mation is somewhat rough due to our limited scan area, 
we expect that the stray field does not exceed 1 Oe, and 
is therefore below the accuracy of our vector magnet with 
a resolution of 1 Oe. 

In order to study the depinning of vortices in 
ErNi 2 B 2 C, we manipulate the separate, stray field in¬ 
duced vortices via the tip-sample interaction. In Figs. 
2(a) and 2(b) two 15 fim x 15 (i m MFM images are 
shown at T = 1.5 K and lift heights ( LH ) of 200 nm 
and 300 nm. The scan area corresponds to the one in 
Fig. 1(a). While at LH = 200 nm a manipulation of the 
lower vortices along the twin boundary is observed, an 
increase to LH = 300 nm leads to an image of static vor¬ 
tices. For all images the fast scan direction is horizontal, 
while the slow scan direction is vertical, from the bottom 
to the top. Note that the direction of the manipulated 
vortex motion always follows the bright stripe. We recall 
that the magnetic flux direction within the twin bound¬ 
ary is anti-parallel to that of the stray field vortices, as 
evidenced by the opposing contrast that correspond to 


opposing magnetic forces acting on the tip. It is therefore 
energetically unfavorable for a separate vortex to move 
into the stripe. On the other hand, the upper right vor¬ 
tex caught between two domain boundaries can move nei¬ 
ther up nor down and hence experiences a much stronger 
pinning. To estimate the pinning force of vortices as a 
function of temperature, we consider the magnetic force 
between tip and vortex, given by F tv = • ^ z+ x(T)] 2 ^ 

Here, m = (3.7 d= 0.2) nAm is the magnetic moment per 
unit length of the tip, <f>o = h/2e is the magnetic flux 
quantum of a single vortex, and z is the distance between 
tip and sample surface around which the vortex manip¬ 
ulation sets in. m has been estimated via monopole- 
monopole approximation. 27 ! Considering that at T = 1.5 
K the onset of the vortex manipulation is at a tip-sample 
distance of 200 nm - 300 nm, we can get a rough estima¬ 
tion of the pinning force Fp for vortices in the range of 
5 - 8 pN. 

When increasing the temperature from 1.5 to 4.5 K and 
thereby moving from the weak ferromagnetic phase into 
the antiferromagnetic phase, the magnetic stripes vanish 
while the twin boundaries remain [see Figs. 2(c) - 2(f)]. 
Hence, the background in the MFM images appears to 
be homogeneous. In contrast, stray field induced vortices 
remain as expected. In this intermediate temperature 
regime the pinning force is strongly reduced compared 
to the weak ferromagnetic phase. This is in good ac¬ 
cordance with previous magnetization and transport ex¬ 
periments, where a threefold decrease in the relative Fp 
was reported when heating the sample up through Twfm 
[see solid blue line in Fig. 2(g)] PS Still, the direction of 
vortex manipulation remains along the twin boundaries, 
which continue to carry a reduced ferromagnetic com¬ 
ponent along the c axis and anti-parallel to the vortex 
flux. This agrees well with the previous Bitter decora¬ 
tion experiment, where field-induced vortex stripes occur 
along [110] and [110]. 14 Their work reports parallel vortex 
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tip lift height 


FIG. 2. (Color online) (a) and (b): 
MFM images taken at T = 1.5 K 
with 200 nm and 300 nm tip lift 
height (. LH ). (c) - (f): MFM images 
taken at T = 4.5 K with 300 nm, 
400 nm, 500 nm, and 600 nm tip lift 
height, (g) The pinning force of the 
vortices obtained from MFM mea¬ 
surements as a function of temper¬ 
ature (red circles) compared to the 
pinning force estimated from mea¬ 
surements of the critical current I c 
(solid blue line) 


stripes with 1-2 /im spacing, which is also confirmed by 
our observations of parallel vortex lines aligned on both 
sides of the twin boundaries [see arrows in Fig. 2(c)]. We 
manage to manipulate the vortex position at 4.5 K up to 
a tip-sample distance of 500 nm. Only at LH > 600 nm, 
a static vortex image is obtained, leading to a pinning 
force of 2.0 - 2.5 pN. 

Moving from the magnetically ordered (T < T/v) into 
the paramagnetic superconducting regime (T > T/v), we 
study the vortex dynamics as a function of temperature 
at a constant LH of 700 nm (see Fig. 3). It is important 
to note that the vortex configuration changes strongly 
around T/v = 6 K [see Figs. 3(a) and 3(b)]. A detailed 
investigation of the critical current revealed a dynamic 
domain formation and annihilation in a narrow temper¬ 
ature regime around T/v, which leads to a considerable 
increase in antiferromagnetic antiphase domain bound¬ 
aries on which vortices strongly get pinned. 29 ^ These do¬ 
main dynamics are evidenced by our observation of vor¬ 
tex rearrangements at T = 6 K. Above T/v, the struc¬ 
ture relaxes from orthorhombic to tetragonal, leading to 
a disappearance of the twin boundaries, and hence to a 
random distribution of vortices. This is contrasted by the 
Bitter decoration study,®’ where a clear hexagonal vor¬ 
tex lattice forms for T > T/v in an applied field of 25 Oe. 
Since our stray field is too weak to induce a considerable 
vortex density, local impurities may dominate the vortex 
arrangement instead. The increasing thermal energy of 
the vortices allows a manipulation with a LH of 700 nm 
at 7 K and higher, as it leads roughly to a twofold re¬ 
duction in the pinning force compared to T = 4.5 K. We 
estimate a pinning force at 7 K of around 1.0 - 1.5 pN. 
The temperature dependence of the pinning force Fp is 
plotted in Fig. 2(g). It mirrors the critical current be¬ 
havior, which is closely related to Fp as both exhibit 
roughly a twofold decrease between 4.5 and 7 K. More¬ 
over, our Fp data is consistent with the observation of 
a finite J c at T > T/v in Ref. which was previously 
conjectured to vanish above the Neel temperature! 30 At 


7 and 8 K in the absence of twin boundaries, there is no 
more preferred vortex manipulation direction. Instead, 
vortices are being dragged equally to the left and right 
while a general upwards trend in motion is observed, thus 
following the slow scan direction. Eventually, thermally 
enhanced vortex dynamics overcome the pinning poten¬ 
tial and ErNi 2 B 2 C transits into a vortex liquid phase 
between 8 K and 9 K, i.e. close to T c . 

Our MFM images also allow us to comment on the 
issue of the coexistence mechanism in ErNi 2 B 2 C. It is 
believed that magnetism and superconductivity can co¬ 
exist under two possible scenarios: (i) a spontaneous 
formation of vortices 31 32 or (it) a variation of magnetic 
moments on a length scale smaller than the penetration 
depth A. 31 33 We clearly detect a modulation of the mag¬ 
netic background over several micrometers. However, ad¬ 
ditional small-scale magnetic variations are not observed 
in our MFM images down to our resolution limit of 20 
nm. Hence we do not find direct evidence for scenario 
(ii) from our data, although we cannot exclude possible 
magnetic variations on a smaller length scale than the in¬ 
strumental limit. In addition, a phase separation between 
purely ferromagnetic and superconducting regions is not 
evident as the magnetic stripes still exhibit a Meissner 
force. On the other hand, scenario (i) is possible if the 
condition H cl < B int < H c2 is fulfilled! ^ For ErNi 2 B 2 C, 
the critical fields H c \ and H c2 are around 500 Oe and 12 
kOe (below Twfm), respectively! 11 As B int = 500 — 700 
G j 11 1 12 we have a system at hand in which the sponta¬ 
neous formation of vortices could very well be realized at 
T < Twfm • We can estimate the number of vortices n 
expected to appear in the presence of the intrinsic mag¬ 
netic field. For a 15 x 15 /im 2 image as shown in Fig. 
1(a) our estimation yields n = Bi £ t Q A ~ 7600 vortices. 
It is tempting to reinterpret the magnetic stripes below 
Twfm as regions of ultra-high vortex density, where sin¬ 
gle vortices cannot be resolved separately anymore. Con¬ 
sidering that the twin boundaries carry a ferromagnetic 
component aligned along the c axisp 0 and that the flux 
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FIG. 3. (Color online) (a) A 15 pm x 15 pm MFM image obtained at T = 5 K, (b) - (e) 20 /im x 20 pm MFM images taken 
at T = 6 K, 7 K, 8 K, and 9 K. The lift height for all images was constant at 700 nm. 


lines of the vortices induced by the weak intrinsic fer¬ 
romagnetism run along the c axis as well, an accumu¬ 
lation of vortices in the vicinity of the twin boundaries 
could be energetically favorable. In this case the locally 
probed A values in the striped region [see black / green 
dot, Fig. 1(a)] do not reflect the true penetration depth, 
as each force gradient curve has a contribution of both the 
Meissner force from the superconducting region as well 
as the magnetic force from the flux quanta. In a previ¬ 
ous Bitter decoration study vortex stripes pinned along 
the [100] direction have been observed at temperatures 
slightly below Twfm and in applied magnetic fields of 
15 - 20 Oe, which rearrange into a hexagonal vortex lat¬ 
tice above Twfm They were interpreted as a pinning 
of external field-induced vortices along ferromagnetic do¬ 
mains, that vanish above Twfm [see Figs. 2(b) and 2(c)]. 
A spontaneous flux lattice, on the other hand, was not 
observed. In our experiments, carried out in negligible 
external magnetic fields, the stripes vanish completely 
above Twfm • This is in stark con trast t o several previ¬ 
ous Bitter decoration experiment^ 14 20 21 and further sup¬ 
port for our scenario of the spontaneous vortex formation 
induced by the weak intrinsic ferromagnetic field. Note 
that also a slight rearrangement of stray field induced 
vortices can be found upon crossing Twfm , which could 
result from an interaction with departing spontaneous 
vortices. In order to unambiguously prove the existence 
of spontaneous vortex stripes in ErNi 2 B 2 C, however, fur¬ 
ther experiments with a higher resolution or within a 
lower vortex density region (e.g., by applying an oppos¬ 
ing magnetic field) are required in a narrow temperature 
regime around the weak ferromagnetic transition. 


IV. SUMMARY 

In summary, we have presented an MFM study of the 
magnetic superconductor ErNi 2 B 2 C down to sub-Kelvin 
temperatures. We introduced a straightforward method 
for locally determining the vortex pinning force Fp, as 
well as demonstrated a spatially resolved measurement 
of the penetration depth A. Our measurements yield 
A(0) ~ 190 nm. Below Twfm , A varies locally by more 
than 20% on a micrometer scale, while following closely 
changes in the intrinsic magnetic field. The pinning force 
of stray field induced vortices shows a strong tempera¬ 
ture dependence that is related to the existence of var¬ 
ious magnetic phases. We estimate the pinning forces 
to be around 1-8 pN, depending on the temperature. 
At the same time we find that Fp varies locally, and 
easy vortex manipulation directions exist in the magnetic 
phases, but are absent above T/v- Finally, a micrometer 
scale magnetic background detected at T=500 mK that 
coincides with the twin domain boundaries and vanishes 
above Twfm supports the scenario of spontaneous vor¬ 
tex stripe formation in ErNi 2 B 2 C. 
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